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2 ■ Abstract 

The spectra of light -light and heavy-light mesons are described by spinless Salpeter 
equation and Dirac equation respectively, which predict linear dependence of the meson 
mass squared M^ on angular momentum J and number of radial nodes n. Both spectra 
are computed by the WKB method and shown to agree with exact numerical data 
^ . within few percent even for the lowest levels. The drawback of Salpeter and Dirac 

equation is that (inverse) Regge slopes do not coincide with the string ones, 27rcr and 
cn ' vro" respectively, because the string dynamics is not taken into account properly. The 

^^ . lacking string rotation is introduced via effective Hamiltonian derived from QCD which 

Qs \ generates linear Regge trajectories for light mesons with the correct string slope. 

Ph! 1 Introduction 

O . QCD is believed to be the fundamental theory of strong interactions and the meson spec- 

troscopy is to be derived from QCD. The spectrum of mesons has been treated in a sequence 
of models [|l], H which may be called QCD motivated, but still not directly derived from the 
/\ ■ QCD Lagrangian. 

c^ ■ The problem of the celebrated Regge behaviour of the hadron spectra has been discussed 

in literature not once (see e.g. |0] and references herein) but still attracts considerable atten- 
tion. The light-light meson spectrum obtained so far and reasonably describing experiment 
can be written in the form 

Mfi{n, J) = (c^'^n + cfj + AM^^ ^ ^^2)^ (1) 

where n is the radial quantum number, J being the total angular momentum, AM^ contains 
the perimeter (self-energy) mass correction as well as corrections to the first two terms, while 
AM^ takes into account spin splittings. 

For heavy-light mesons a similar relation holds true with the subscript // changed for hi 
in all coefficients: 
M,^,(n,J) = ( ci^')n + cf^J + AMI + ^^')- (2) 
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From general physical considerations one expects the spectra (^ and (0) to follow from 
a string-like picture of confinement which predicts the (inverse) Regge slope 

cf = 27r(7, cf ) = Tia. (3) 

Additional daughter Regge trajectories are given by vibrational excitations missing in 
(H]) and (H), which are due to hybrid excitations, i.e. constituent gluons attached to the 
fundamental string 0] . In what follows we are interested only in radial and orbital excitations 
of the string. 

The string slope @) is an important criterion and check for any QCD inspired model 
since it requires the correct account for the rotation of the string, which is not present in 
the potential models considered so far. For example, relativistic spinless Salpeter equation 
with confinement reduced to the linearly rising potential between quarks yields 

Cj {Salpeter) = 8o", Cj {Salpeter) = Aa (4) 

that is about 25% larger than (|^), whereas the one-body Dirac equation with the linearly 
rising potential leads to 

cf\Dirac) = 4:a, (5) 

if the potential is added to the energy term (vector confinement[]) , and 

cf\Dirac) = 2a, (6) 

for the potential added to the mass term (scalar confinement). Both results lead to consid- 
erable discrepancies with (|^) and, as will be shown later, this happens because the rotation 
of the string, and hence momentum dependence of the effective potential, is not taken into 
account. 

It was found few years ago |^ that starting from the area law for Wilson loops one arrives 
at the relativistic Hamiltonian for the spinless quark and antiquark which possesses two 
different regimes: potential regime for small angular momenta L and any n, and string-like 
one for large L and fixed n. In the latter case the dominant term in the Hamiltonian indeed 
describes the rotating QCD string, so that the string Regge slope (|^) is readily reproduced. 
Similar results were obtained independently by numerical analysis of the spinless quark- 
antiquark system 0. 

In the present paper we concentrate on the quasiclassical approach to mesons, as the 
WKB method allows to obtain analytic formulae for the meson spectra of surprisingly high 
accuracy thus giving evidence for the quasiclassical dynamics of confined quarks in the meson. 

Therefore our first task will be to check the accuracy of the WKB approximation for 
those cases where exact solutions are feasible: spinless Salpeter equation for light-light 
mesons (potential regime of the general QCD string formalism 0) and the Dirac equation 
for linear confining potential for the case of heavy-light system. We argue that the accuracy 
of WKB results is very good even for lowest states. However the slopes in both cases are 
incorrect, as in (^ and (|) respectively. 

At this point we come to the main purpose of this study — to include the proper string 
dynamics, whereby abandoning the notion of local potential and introducing a new entity. 



^Here we leave aside the well-known problem of the Klein paradox revealing itself in case of vector 
confinement. 



the QCD string, the effect which can not be recasted in terms of local potential. We use the 
Hamiltonian derived in [^ and calculate the quasiclassical spectrum of light mesons. The 
results represent celebrated straight-line Regge trajectories even for low-lying states with the 
slope very close to the expected string slope (||). 

In conclusion we demonstrate how other effects (spin and colour Coulomb interaction) 
can be included in the same Hamiltonian to make a direct comparison with experiment. 

2 Meson spectrum and quasiclassical approximation 

We start with the spinless Salpeter equation which describes relativistic quark and antiquark 
of equal masses m with angular momentum / = and spin effects neglected (see for the 
derivation of this equation from the general meson Green function in QCD). 

(2^p2 + ^2 + ^^)^^ = jvff )^„ (7) 

The Bohr-Sommerfeld condition looks like 

J^ Mr)dr = n[n + ^y n = 0,l,2,..., r+ = ^^^^^^— ^^ (8) 

that yields 

Mf>V(M.?")^-4,n^- 4,n^i,y'"^"''^;^'"^+ig = 4a. (n + 5) (g) 

A similar consideration for the heavy-light system of masses m and M (M -^ cxd) gives 



(Vp^+^+ ar)^„ = Mf '^V-n (10) 

Mf ')y' (^""'0 -'^'- ^''l^— —^ — = 2^^ r + i) • (11) 

Accuracy of WKB approximation (^, (|TT|) can be tested vs exact solutions of the Salpeter 
equations (recently accuracy of WKB approximation was checked for light-light mesons in 
[^). In Table 1 this comparison is given for the light-light system with m = and heavy- 
light one with nig = 0.01 GeV and and Mg = 10 GeV. The mass M^'^') in the latter case 
actually refers to the difference of the total mass of the heavy-light system and the mass of 
the heavy antiquark. 

Summarizing, one can say that spectra (|), ([TT|) (as function of n for / = 0) indeed have 
the form (|I]), (Q) with the corrections at large n in the form 

AM^ = O l^ln^] = O f ^1 . (12) 

The WKB spectrum is linear in n and its accuracy is about 3-4% even for the lowest 
state. 



n 





1 


2 


3 


4 


5 


M^"HWKB) 


1.373 


2.097 


2.629 


3.070 


3.455 


3.802 


M^") (exact) 


1.412 


2.106 


2.634 


3.073 


3.457 


3.803 


Mf'^(WKB) 


0.971 


1.483 


1.859 


2.171 


2.443 


2.688 


Mf') (exact) 


1.014 


1.524 


1.917 


2.246 


2.537 


2.800 



Table 1: Comparison of the WKB and exact spectra for Salpeter equations (0) for massless 
quark and antiquark and (p!0| ) for the quark mass m^ = 0.01 GeV and antiquark mass 
M„- = 10 GeV. In both cases a = 0.2 GeV^ and / = 0. 



We now turn to the case of the Dirac equation with hnear confining potential studied in 
[^. The WKB method for the Dirac equation was thoroughly investigated in [^j and recently 
applied to the case of confining potential [|l^]. Let us briefly recall the results here. 

The Dirac equation with scalar (f/) and vector {V) local potentials has the form 



{dp + (3{m + U) + V)il)n = en-iJJn, 
and the WKB quantization condition is [Q 

r 



\p-\ ] dr 

r- \ pr J 



7T [n -\- 



n 



0,1,2, 



where 



P 



{e-VY---{m + U)\ 



(13) 



(14) 



(15) 



w 



1 
2r 



U' ~V' 



2m + U + e-V' 
1 



J + 



An approximate quasiclassical solution of (|T^) obtained in ||T0[ for the case m = 0,V = 0, 
U = ar is 



2a\2n + j + ^+-^ 



'0.38 + /n^^ 



Txet 



(y\K\ 



+ 




(16) 



The last two terms on the r.h.s. of equation (|T^) are sub-leading for large n and are 
generated by the term ^^ (see (0)). One can see that the (inverse) Regge slope in j in 
(|TBp is equal to 2cr coinciding with the exact result (|^), but is not of string type. As it 
was expected a j-independent scalar potential does not describe the physical phenomenon 
of rotating string. 

Still the accuracy of the WKB approximation is impressing. In Table 2 one can see the 
comparison of exact eigenvalues computed in |^ with quasiclassical ones and with those 



obtained from (16). The discrepancy is less then 1% even for the lowest state and it is much 
better for higher states. 
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J 
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Mexact 


MwKB 


^eq.(R) 





2 


3/2 


2 


1.209 


1.209 


1.208 





1 


1/2 


1 


1.024 


1.025 


1.025 








1/2 


-1 


0.725 


0.726 


0.733 





1 


3/2 


-2 


0.957 


0.960 


0.966 





2 


5/2 


-3 


1.149 


1.149 


1.155 



Table 2: Comparison of the exact spectrum for Dirac equation ([T3|) with WKB spectrum and 
approximate one given by equation (|1^) for m = 0, V^ = and U = ar with a = 0.2 GeV^. 



3 Rotating string in the spinless quark Hamiltonian 

Let us turn back to the spinless Salpeter equation for the light -light meson and take the 
non-zero angular momentum into account. As it was shown above the Salpeter equation with 
local /-independent potential leads to the incorrect Regge slope (Q), and therefore this case 
requires a special treatment. One needs a Hamiltonian taking into account dynamical degrees 
of freedom of the string, e.g. in the form of time derivatives of string coordinates. This was 
done explicitly in [Q , where it was shown that starting from the QCD Lagrangian and writing 
the gauge invariant qq Green function for confined spinless quarks in the Feynman-Schwinger 
representation, one can arrive at the Lagrange function of the system in the well-known form 



L(t) 



-mi 



m2i 



a 



df3 



ww'Y — w'^w''^. 



(17) 



where r denotes the proper time of the system, the first two terms stand for quarks, whereas 
the last one describes the minimal string with tension a developed between the constituents; 
w^{t,(3) being the string coordinate. Adopting the straight-line anzatz for the minimal 
string, i.e. w^{t,P) = /3xi^ + {1 — P)x2fj., synchronizing the quarks proper times, Xio = X20 = 
T = Uah and introducing auxiliary fields to get rid of the square roots (see e.g. |]ll|) one can 
obtain the following Hamiltonian in the centre of mass frame (we consider the case of equal 
masses m) 



H 



pI + rri^ 



Kv 



+ M^) + 



a^r^ 






^ 2 Jo v{l3) ' Jo 2 

where the two auxiliary positive functions fi{T) and z/(/3, r) = i^(/9) are to be varied and to be 
found from the minimum of H yielding quark energy and string energy density respectively. 
A more detailed analysis of the role played by auxiliary fields can be found in e.g. [T2\. 
Note that Hamiltonian (II8D has a form of the sum of "kinetic" and "potential'' 



terms 

only due to auxiliary fields /i and u. If one gets rid of them by substituting their extremal 
values, the resulting Hamiltonian possesses a very complicated form which makes its analysis 
and quantization hardly possible. 

The centrifugal potential in Hamiltonian (|THp is of special interest to us and, most of 
all, the second term in the denominator. It is this term that describes extra inertia due to 
the string connecting the quarks. Neglecting this term and taking extrema in the auxiliary 
fields one easily arrives at the ordinary Salpeter Hamiltonian with linearly rising potential. 



whereas account for this extra term describes the proper string rotation and brings the slope 

of the Regge trajectory into correct form (|]). In the nonrelativistic expansion of Hamiltonian 

18) this term yields the so-called string correction to the leading confining potential ar |^ 



Aff, = ""^ 



Qrn?r' 



the part of the interaction which explicitly depends on the angular momentum. 

Hamiltonian (|l^) assumes especially simple form in the case of zero angular momentum 
and after excluding the auxiliary fields produces Salpeter equation (0). 

Variation of (0) over z/(/3) gives the stationary energy distribution along the string with 
/? (0 < /? < 1) being the coordinate along the string. Thus one obtains 

i/o(/5) = '"'' (19) 

^l-V(/3-l)2' 

where y is to be found from the transcendental equation 

-— = -—{arcsin y - yJl-y^) + — , (20) 

UJ.Z /^yZ V uj. 

andL2 = /(/ + l). 

Note that the maximal possible value of y, y = 1, yields the energy distribution z/q ^^(/5) 
corresponding to the free open string (string without quarks at the ends) P, ||. 

In the general case inserting the extremal function vq{j3) one obtains from ( |18|) 



H = "'' h /x(r) + '^arcsin y + \i{r)y^ (21) 

/^(^) y 

with y = y{L, r, /x) defined by equation (|^). Unfortunately no rigorous analytic calculations 
are possible anymore, so one has to rely upon numerical calculations. But let us first perform 
some analysis of Hamiltonian ( pT]) . 

Neglecting fi in ( PP] ) and fiy"^ in ( pi| ) (which is justified for large L and ar, so that — ^ 1) 
and varying over /i in (pT]) one obtains 



ar 



Has = '^yPr + rn? -\ arcsin y, (22) 

so that the second term on the r.h.s can be viewed as an effective potential, and we would 
like to emphasize that this potential is non-trivially /-dependent. 

In the general case one has a /^-dependent Hamiltonian (|T8|) with the "potential" U{fj,, r), 

U(fi, r) = — arcsin y + yuy^. (23) 

y 

A simplifying approximation can be used at this step, namely the standard WKB proce- 
dure can be applied to the Hamiltonian 

9 o 

H = ^^-^^ + fio + U{fio,r), (24) 



which shghtly differs from the exact Hamiltonian ( [2T| ) as it treats fiQ as a variational param- 
eter not depending on r. We find eigenvalues M{fiQ,L,n) and minimize them with respect 
to fiQ to obtain the spectrum M{fiQ{L,n),L,n), where nl^L^n) being the extremal value of 
/io- 

To check the accuracy of such a procedure for the eigenvalues two Hamiltonians were 
considered: 

Hi =2^pl + m^ + ar, (25) 



H, 



pI + m^ 

1^0 



+ fJ'O 



or. 



/io varied, 



(26) 



where Hi is obtained from H2 in the limit when /xq — ^ /u(t)- 

The results are listed in Table 3. One can see that the accuracy of variational procedure 
(pBD is better than 5% and it is reasonable even for m tending to zero. 



n 





1 


2 


3 


4 


5 


MU) 


1.475 


2.254 


2.825 


3.299 


3.713 


4.085 


Mf) 


1.412 


2.106 


2.634 


3.073 


3.457 


3.803 



Table 3: Comparison of the WKB spectrum of Hamiltonian (^61) M^^^ with the exact spec- 
trum of Hamiltonian ( P^D M^'> for m = and a = 0.2 GeV'^. 



As a next step we use the standard WKB method to find the spectrum of Hamiltonian 
41). To this end we write the Bohr-Sommerfeld condition as 



''+ / 1 

Pr{r)dr = 71 in + - 

T- V 2 



(27) 



with 



Pr{r) = ^/io(M-/io- t/(/io,r)) 



m^ 



:28i 



The eigenvalues M{fiQ,L,n) were found numerically from (^) and the minimization 
procedure was then used with respect to /iQ. Results for Mni are given in Table 4 and 
depicted in Fig.l demonstrating very nearly straight lines with approximately string slope 
(27ro")~^ in / and as twice as smaller slope in n. 



n I 


1 


2 


3 


4 


5 





1.865 


2.200 


2.481 


2.729 


2.956 


1 


2.562 


2.832 


3.068 


3.281 


3.480 


2 


3.091 


3.329 


3.540 


3.733 


3.913 


3 


3.535 


3.753 


3.947 


4.125 


4.290 


4 


3.925 


4.128 


4.309 


4.476 


4.629 


5 


4.278 


4.469 


4.638 


4.797 


4.939 



Table 4: Quasiclassical spectrum of Hamiltonian (0) for m = and a = 0.2 GeV^. 

Let us give a little comment concerning effective potential U{^o,r). Its behaviour at 
large and small distances can be extracted analytically from Hamiltonian (|T8|) and coincides 
with that of the Salpeter: the centrifugal barrier at small r and linear growth at large r. 



Meanwhile in the region of intermediate values of r this potential differs from what one 
would have in the Salpeter equation and it is just this region which is important to obtain 
the correct Regge slope. The form of the effective potential is depicted in Fig. 2 for several 
different angular momenta /. In case of / = the effective potential equals to ar for all 
values of r. 



4 Conclusion 

We have shown that the proper account of the string dynamics leads to practically linear 
Regge trajectories, shown in Fig.l, with the slope numerically close to the conventional 
(27ro")~^. The exact form of the effective potential incorporating the string rotation as well 
as the quark radial motion was found numerically and shown in Fig.2. 

To make contact with experimental data on meson masses one should specify correc- 
tions AMp, AM^ in (|l]), or in the case of the Hamiltonian formalism, one should add to 
Hamiltonian ( PT| ) the colour Coulomb term Vq and spin-dependent interaction. 

Treating the latter as perturbation one finds, e.g. for p meson, a negative shift of the mass 
due to Vc of about 160 MeV and positive correction of about 40 MeV due to hyperfine term 
AHss- Taking this into account one obtains for p meson (/ = 0) the mass about 1.6 GeV. 
It is clear from Fig.l that these corrections practically do not violate the linearity of Regge 
trajectories as p meson lies on the continuation of the leading theoretical trajectory in I 
(see dashed line attached to the trajectory with n = 0). In this way starting from QCD 
and making one assumption of the area law for the Wilson loop we obtain linear Regge 
trajectories for light quark mesons with the string slope. 

In this discussion quark spin effects have been taken into account perturbatively, which 
is a reasonable approximation for the p trajectory, but unacceptable for pions and kaons, 
since for the latter one needs the full implementation of the chiral dynamics. 

The progress in this direction was achieved in recent papers of one of the authors (Yu.S.) 
p!0[] , where an effective Dirac equation for the quark moving in the field of an infinitely 
heavy antiquark source was derived, and it was shown that solutions display the proper- 
ties of confinement and chiral symmetry breaking. The nonrelativistic limit of the resulting 



interaction lead to the conventional result for the confining term and to the spin-orbit in- 



teraction in agreement with the standard Eichten-Feinberg-Gromes results [|T3]. Therefore 
pionic trajectories should be considered in this new formalism. 

There is yet another question unanswered by our paper (and to our knowledge by all 
other existing papers) — the intercept of Regge trajectories Lq = L (M^ = 0). Theoretical 
intercept for the leading trajectory in j (see Fig.l and the caption to it) is around -0.5, 
whereas it is -1-0.5 for the experimental p trajectory also shown in Fig.l. The customary 
way in the potential models is to add to the Hamiltonian a large negative constant |Co| ~ 
1 GeV to reproduce the intercept, but this would obviously violate the linearity of Regge 
trajectories. Therefore one expects that QCD provides a negative constant AM^ in (|1|) but 
not in Hamiltonian (|l^). 

The authors are grateful to A.M.Badalian, A.B.Kaidalov, Yu.S.Kalashnikova and V.S.Po- 
pov for useful discussions. Financial support of RFFI through the grants 97-02-16404, 97- 
02-17491 and 96-15-96740 is gratefully acknowledged. 
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Figure 1: Quasiclassical spectrum of Hamiltonian ([1^) for m = and cr = 0.2 GeV^. The 
leading experimental Regge trajectory in angular momentum / is given in the upper plot for 
comparison. Note that as soon as / = j — 1 for the given experimental trajectory, then it lies 
somewhat higher when plotted in / than when plotted in j. Theoretical prediction for the 
p-meson mass, M^ ^ 2.5 GeV^, with colour Coulomb interaction and spin effects included, 
is shown not to violate the straight-line behaviour of the leading theoretical trajectory. 
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Figure 2: Effective potential incorporating the string rotation as well as the quark; radial 
motion for a = 0.2 GeV^. Dashed lines give matching with the centrifugal barrier at small 
r, whereas for all / > the effective potentials have the same asymptotics ar at large r. For 
/ = the effective potential coincides with ar for all values of r. 
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